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ABSTRACT 

We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the 
stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the 
south-west semi-major axis. These data enable the construction of colour-magnitude diagrams 
reaching the oldest main-sequence turn-offs (^13 Gyr) which, when combined with another 
field at 25 kpc from our previous work, we use to derive the first precision constraints on 
the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit 
temporal as well as field-to-field variations, but are generally always within a factor of two 
of their time average. There is no evidence of inside-out growth over the radial range probed. 
We find a median age of ^7.5 Gyr, indicating that roughly half of the stellar mass in the M31 
outer disc was formed before z 1. We also find that the age-metallicity relations (AMRs) 
are smoothly increasing from [Fe/H]~-0.4 to solar metallicity between 10 and 3 Gyr ago, 
contrary to the flat AMR of the Milky Way disc at a similar number of scale lengths. Our 
findings provide insight on the roles of stellar feedback and radial migration in the formation 
and evolution of large disc galaxies. 

Key words: galaxies: individual: M31 - Local Group - galaxies: formation - galaxies: evo¬ 
lution - galaxies: structure - galaxies: stellar content 


1 INTRODUCTION 

Disc galaxies account for a sizeable fraction of the stellar mass in 
the Universe, yet the details of their formation and evolution are 
still greatly debated. While much progress has been made in simu- 
lating realistic disc galaxies within a cosmological frame work (e.g. 
lAumer et alj2013l : lMarinacci. Pakmor, & Springel'TOf^ . many is¬ 
sues remain unresolved. Strong stellar feedback at early epochs is 
often included in hydrodynamical simulations in order to suppress 
high redshift star formation, postponing the formation of extended 
discs until relatively recent times (z<l or ~8 Gyr) when the bulk 
of the merging is over fe.g. lStinson et alj|2013h . This scenario nat¬ 
urally results in the inside-out growth of discs, predicting radial 


* Based on observations made with the NASA/ESA Hubble Space Tele¬ 
scope, obtained at the Space Telescope Science Institute, which is operated 
by the Association of Universities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. These observations are associated with pro¬ 
grammes GO-9458 and GO-13018, 
f E-mail: ejb@roe.ac.uk 


age gradients and young mean stellar ages at large radii. Direct 
observations of old and intermediate age stars in outer discs are re- 
quired to shed light on the efficacy and role of e arly feedback (e.g. 
iFerguson & Johns^l200ll : |Williams et alj2013l) . 

Once in place, various processes may rearrange the spatial dis¬ 
tribution of the stellar populations within galactic discs. For exam¬ 
ple, it has been shown that stars can undergo large radial excursions 
within di scs due to scattering off transient spirals and o ther fea¬ 
tures fe.g. lSellwood & Binnew200l : lMinchev et al.l201lh . By effi¬ 
ciently mixing stars of different ages and metallicities over kilopar- 
sec scales, this process could lead to a smearing of the star forma¬ 
tion history (SFH) and age-metallicity relation (AMR) at a given 
radius, therefore complicating the interpretation of the resolved fos¬ 
sil record. Determining the importance of this process in real galaxy 
discs is obviously of paramount importance. 

To explore these issues, we have obtained deep Hubble Space 
Telescope (HST) pointings with the Advanced Camera for Surveys 
(ACS) in order to derive the detailed SFH and AMR at several lo¬ 
cations in the outskirts of M31. The SFHs that we calculated for 
the first 14 fields revealed a global burst of star formation in the 
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Figure 1. Location of our HST/AC S pointings superi mposed onto the 
INTAVFC map of M31’s inner halo iirwin et ^l2005h . showing the dis¬ 
tribution of evolved giant stars around M31. The three fields analyzed in 
this work are shown as open squares, while the open circles represent some 
of the fields studied in Papers I and II. The ellipse has a semimajor axis of 
2° (27 kpc) and represents an inclined disc with i = 77? 5 and position angle 
of38?l. 


M31 outer disc roughly 2 Gyr ago, and allowed us to determine 
the nature and origin of the various subst ructures in the outer disc 
and inner halo l lBemard et al.ll201^120151 hereinafter Papers I and 
II, respectively). In this Letter, we present the SFHs calculated for 
three very deep fields located between 20 and 25 kpc from the cen¬ 
tre of M31 along the south-west semi-major axis. Based on colour- 
magnitude diagrams (CMDs) that reach back to the oldest main- 
sequence turn-offs (MSTO), these data provide the first precision 
view of the spatially-resolved SFH in any disc galaxy and lead to 
a number of interesting conclusions about the assembly history of 
the M31 outer disc. The paper is structured as follows: in Sectionj^ 
we present the observations and the data reduction steps, while the 
CMD-fitting method is briefly described in Section The result¬ 
ing SFHs and our interpretation of the results are discussed in Sec¬ 
tion |4] 


2 OBSERVATIONS AND DATA REDUCTION 

The 20 and 23 kpc fields were observed with the ACS onboard the 
HST (proposal GO-13018, P.I.: A. Ferguson) between 2013 May 
and August. As the goal of the observations was to study the his¬ 
tory of the outer thin disc, the fields were placed very close to 
the sem i-major axis. When combined with the Warp fiel43 anal¬ 
ysed in IPaper ll the three fields probe a range of galactocentric 
distance from 20.1 to 25.4 kpc along the southwest thin disc, cor¬ 
responding to -^3.8 to 4.8 scalelengths in the M31 disc assuming 


Rd = 5.3±0.5 kpc iCourteau et alj2oTll) . In terms of radial scale- 
lengths, we note that the radius of our innermost field is comparable 
to that of the solar neighbourhood Our fields are shown as open 
squares in Fig.[T] overplotted on a map of the surface de nsity of red 
rant b ranch (RGB) stars from the INTAVFC survey of lirwin et al.l 
20051) . The HST field lying on the maj or axis ju st inside the 27 kpc 
ellipse is the Outer Disc field studied in lPaner iL for which we could 
not calculate the SFH because of strong differential reddening from 
dust within M31; the other open circles represent some of the in- 
ner halo and substructure fields for which SFHs were derived in 
IPaper III While the Warp was observed for 10 HST orbits, we ob¬ 
tained 13 orbits per field for each of the new fields in order to com¬ 
pensate for the effect of higher stellar density and background flux. 
A detailed summary of observations including coordinates, fore¬ 
ground color excess, radial distances, and exposure times is given 
in Tables 

The photometry and artificial star tests we re carrie d out fol¬ 
lowing identical methods as those described in IPaper jj The only 
difference here is that we took advantage of the images corrected 
for the charge tran smission inefficiency (‘FLC’) provided by the 
HST pipeline (e.g. lAnderson & Bediiil l2()10l : lUbeda & Anderso^ 
l2012h . The stellar photometry was carried out on the individ¬ 
ual exposures with the standard D AOPHOT/ALLSTAR/ALLFRAME 
suite of programs dStetsonl I 1994 I) . They were then cleaned of 
non-stellar objects by applying cuts on the photometric param¬ 
eters given by ALLFRAME, namely the photometric uncertainty 
(ct $1 0.3) and the sharpness, describing how much broader the 
profile of the object appears compared to the profile of the PSF 
(|SHARP| ^ 0.3). After cleaning, the CMDs contain 353 918 and 
296731 stars, respectively, compared to 185 335 in the Warp with 
the same quality cuts. The final photometry was ca librated to the 
VEGA MAG system following the prescriptions of ISirianni et al.l 
ll2005h and the revised ACS zeropointfli then to absolute magni¬ 
tude by_taking_hito_accounLthe distance to M31 ((m—M)o=24.47; 
e.g. McConnach ie_et_alJ|2005b, and the f oreground reddening from 
the lSchlegel, Finkbeiner & Davisi ( [l99^ extinction maps togethe r 
with the updated coefficients from ISchlaflv & Finkbeineil ( 1201 ll) . 
Note that the CMD-fitting technique used in this paper minimizes 
the impact of the uncertainties in distance, mean reddening, and 
photometric zero-points on the solutions by shifting the observed 
CMD with respect to the artificial CMD; therefore, any small er¬ 
rors in these quantities do not affect the results. 

While great care was taken to place t he fields in areas of low 
dust extinction by using the HI maps from iBraun et alj ( l2009l) . we 
found that the 23 kpc field is affected by mild differenti al redden ing 
due to dust in M31. We used the same method as in IPaner j for 
the Outer Disc field to quantify it. We were not able to correct the 
areas with measureable differential reddening, but found that ~40 
percent of the ACS field was not affected. In the following we only 
use the ~ 120 000 stars that reside in the unaffected areas. 

The resulting CMDs for the two new fields plus the Warp 
are shown in Fig. only 120000 stars per CMD are shown 
to ease the comparison. Selected isochrones and a theoreti- 
cal zero-age horizontal -branch (ZAHB) from the BaSTI library 
dPietrinfemi et al.ll2004h are shown as labeled in the inset, both to 
indicate the location of the MSTOs and as a guide for comparing 


^ This field lies in the southern stellar warp, which bends strongly away 
from the major axis. As a result, its true radius in the disc could be somewhat 
larger than its projected radius. 


^ The Sun is located ~3 .6 scalelengths from the Galactic centre, a ssuming 
Rrl.M W ~ 2.2 kpc (e.g. lBovv & Rixll2013l) and Ro = 8 kpc (e.g. lMall^ 

I2OI3I1 . 

^ http://www.stsci.edu/hst/acs/analysis/zeropoints. 
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Figure 2. CMDs f or the 20 kpc (left), 23 k pc (middle), and Warp (right) fields, where selected isochrones (solid lines) and a ZAHB (long-dashed line) from 
the BaSTI library IPietrinfemi et al.ll2004h are overlaid and labelled in the inset. Only 120000 stars per CMD are shown. The error bars show the mean 
photometric errors as a function of magnitude. The contour levels correspond to [20, 40, 60, 80, 100, 120] X 10® stars mag“®. The 25% completeness limit 
from the aifificial stars tests is shown by the dashed lines. 


Table 1. Summary of Observations 


Field 

20 kpc 

23 kpc 

RA (J2000) 

00:38:05.8 

00:37:23.8 

Decl. (J2000) 

+40:05:35 

+39:54:00 

R^(kpc) 

20.1 

23.3 

E(B-v£ 

0.053 

0.055 

Dates 

2013 May 28-Jun 2 

2013 Jul31-Aug 10 

tF606W,^(s) 

4x1149+8x1245=14 556 

4x1145+8x1240=14 500 

tF814w'^ (s) 

4x1149+10x1245=17 046 

4x1145+10x1242=17 000 


“ Galactocentric distance. 

^ Values from ISchlaflv & Finkbeineil i201 ih . 
Individual and total exposure times. 


the CMDs. The 25% completeness limit from the artificial star tests 
is represented as dashed lines. It shows that while the 23 kpc and 
Warp field have the same photometric depth, the 20 kpc field is 
slightly shallower due to the higher stellar density and sky back¬ 
ground at the time of the observations. In terms of stellar pop¬ 
ulations, the main differences between these fields are the grad¬ 
ually less populated main-sequence (at Meoe — Msn ~ 0 and 
M 814 < 2) with increasing radius, and the RGB of the Warp above 
M 814 — 1 being on average bluer by 0.05 mag than in the other 
fields. The Warp also contains a densely populated plume of stars at 
Meoe — Msia ^ 0.4 and M 814 1.5 that corresponds to a ~2 Gyr 

old stellar population (see below), and which is not as prominent in 
the other fields. 


3 STAR FORMATION HISTORY CALCULATION 

The SFHs were calculated using the technique of synthetic CMD 
fitting, which provides a detailed and robust quantitative estimate 
of the star formation rate (SFR) and metallicity ([Fe/H]) at a given 
time. This involved fitting the observed data with synthetic CMDs 
to extract the linear combination of simple stellar populations (SSP) 


- i.e. each with small ranges of age (^1 Gyr) and metallicity 
(<0.25 dex) - which provides the best fit; the amplitudes of the 
required SSPs give the rates of star formation as a function of age 
and metallicity. 

The methodology is identical to that used in Papers I and II, 
although we are now using a new algorithm developed in Python 
by one of us (EJB). The only notable difference is that the good- 
ness of fit is m easured through a Poisson equivalent of p (see 
iDolnhinl l2002h rather than the modified statistics of iMighelil 
1 1999h . Extensive testing has shown that it produces results that are 
virtually indistinguishible from those obtained with the algorithm 
used in our previous papers (lAC-POP/MINNIAC; e.g. lHidalgo et aD 
I 2 OI ih . This can be verified by comparing the SFH of the Warp that 
w as calcu lated using the new algorithm (see Eig. with Fig. 6 
in IPaner E The advantages of the new approach are twofold: the 
algorithm is at least an order of magnitude faster - a necessary effi¬ 
ciency in view of the ever-increasing size of star catalogues - and it 
is significantly easier to use. We therefore refer the interested reader 
to Papers I and II for a detailed description of the assumptions and 
fitting procedure. 


4 RESULTS AND DISCUSSION 

The resulting SFHs are presented in Fig. where the top and 
bottom panels show the evolution of the SFR and of the metal¬ 
licity, respectively. The redshift scale was constructed as suming 
the WMAP7 cosmological parameters jjarosik et alj|201 ih . In all 
fields, we find that star formation began early on (z > 2) and oc¬ 
curred more or less continuously across the history of the disc. This 
indicates that M31 was already a large disc galaxy at high redshift 
and could be a direct descendent o f the large disc systems iden¬ 
tified in distant redshift studies (e.g. iFdrster Schreiber et al.ll2009l ; 
Ivan Dokkum et alj2013l) . The SFHs exhibit complex temporal vari¬ 
ations; aside from the strong decline over the last 2 Gyr, these vari¬ 
ations are not clearly correlated from field to field. Nonetheless, the 
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Figure 3. Evolution of the SFR (top) and metallicity (bottom) as a function of time for the three fields, displayed on the same absolute scale. The dashed lines 
in the top panels show the average SFR over the age of the disc. In the bottom panels, the grayscale is proportional to the stellar mass formed in each bin. The 
filled circles and error bars show the mean metallicity and standard deviation in the age bins representing at least 1 percent of the total mass of stars formed; 
the solid blue line serves to guide the eye. 


SFR in a given field remains within a factor of 2 of its time av¬ 
erage, indicating that very strong bursts of star formation have not 
dominated the mass build-up in the stellar disc. As expected from 
the decreasing stellar density with increasing galactocentric radius, 
the average SFR decreases by a factor ~3.3 over the radial range 
spanned by our data (roughly one scalelength). 

A particular feature that stands out in the SFH of the Warp is 
the burst of star formation -^2 Gyr ago, associated with the onse t of 
a decrease in global metallicity, that we first detected in lPaper j and 
later found to be widespread in more distant (R> 25 kpc) fields in 
M31 jPaper Ilh . Interestingly, while the decrease in [Fe/H] is clearly 
measured in the two new fields presented here, we do not detect any 
clear enhancement in the SFR at that epoch. The 20 kpc field does 
show an enhanced period of star formation activity but this is not 
as pronounced as the Warp burst, and occurs at earlier epochs (3-5 
Gyr), so it is not clearly related. We thus conclude that either the 2 
Gyr burst was confined to the very outer disc, or that it had much 
greater significanc e there. These findings appear to be at odds with 
SFHs calculated bv lWilliams et al.l ( 1201 Sl) in the northern inner disc 
of M31, which show a ubiquitous strong burst 2-4 Gyr ago in fields 
ranging from 3-20 kpc. The reason for this discrepancy is presently 
not clear, although we note that their result is largely based on mod¬ 
elling evolved stars (red clump and RGB) as age indicators, which 
is not as robust as modelling the main sequence turn-offs in our 
much deeper data. 

To facilitate comparison. Fig. |4] combines the SFHs of our 
three fields, where the top and bottom panels show respectively 
the cumulative mass fraction and the AMR as a function of time. 
The top panel reaffirms that even though the details of the stellar 
mass growth vary between the fields, all three show significant star 
formation at « > 1. Combining all three fields, the median stel¬ 
lar age in the outer disc fields is ~ 7.5 Gyr. To test the reliability 
of these old ages, we ran constrained SFH calculations where only 
model populations younger than a given age (8 Gyr, 10 Gyr) were 


allowed. While this affected the ability of the algorithm to fit the 
observed CMD, the median age of the populations did not change. 
With the caveat that our fields span only one scalelength in radius, 
there is no signal of inside-out growth in our data. The epoch at 
which 50% of the present-day stellar mass was in place is 6.9, 8.2, 
and 8.6 Gyr for the 20 kpc, 23 kpc, and Warp fields, respectively. 
The mild trend of older mean age for increasing radius in M31 is in 
cont rast to some recen t integrated light analyses of local galaxies 
(e.g. lPerez et al.|[2oT3h . 

The bottom panel of Fig.|4]shows that, for the past 9-10 Gyr, 
the AMRs of the three fields are very similar within the uncertain¬ 
ties. All start from a moderate level of pre-enrichment and all ex¬ 
hibit a clear decline in the last 3 Gyr. However, we find that the 
older populations of the Warp are considerably more metal-poor 
than in the 20 and 23 kpc fields, consistent with the slightly bluer 
RGB in that field (see Fig.j^. This variation is in keeping with the 
mild abundanc e gradient that has b een measured in this part of the 
M31 disc (e.g. lKwitter et al.ll2012h . 

However, the most striking feature of these AMRs is the rel¬ 
atively tight, smoothly increasing metallicity between ~10 and 
3 Gyr ago. This is in stark contrast with the flat AMR of the so¬ 
lar neighbourhood over the same age range (e.g. lHavwoodl[200^ : 
ICasagrande et al.l 1201 ll) . This comparison is relevant since our 
20 kpc field in M31 is located at roughly the same number of ra¬ 
dial scalelengths as the Sun in the Milky Way. We comp are the 
Geneva-Copenhagen Survey (GCS; l&sagrande et alj20nE AMeQ 
with ours in the bottom panel of Fig.|4j it shows that their analysis 
revealed no metallicity increase in the past 10 Gyr, compared to the 
enhancement of about 0.4 dex that we have measured in the disc 
of M31. We note that more recent works on the A MR of the solar 
neighbourhood based on high resolution spectra dHavwood et al.l 


^ Where the age is their “BaSTI Max likelihood age”. 
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Figure 4. Comparison of the SFHs of the three fields, showing the evolution 
of the cumulative mass fraction (top) and metallicity (bottom), as labeled in 
the inset. The integrated stellar mass density in each field at the present-day 
is 3.9, 2.1 and 1.2 X 10 ®Mq kpc“^ for the 20 kpc, 23 kpc, and Warp 
fields, respectively. In the bottom panel, the contours con'espond to the l-cr 
envelope for the age bins representing at least 1 percent of the total mass of 
stars formed, shown as blue error bars in Fig. [3] The black filled circles with 
error bars show the AMR of the solar neighbourhood (mean metallicity and 
standard deviation) from a reanalysis of the Geneva-Copenhagen Survey 
iCasagrande et al.l201 ih for comparison. 


I 2 OI 3 I : iBergemann et al.ll20l^ do find a clear metallicity increase 
before about 8 Gyr ago, but a similarly fiat AMR since 2 : ~ 1. 

The difference between the AMRs of the M31 and the Milky 
Way discs may indicate different formation and/or subsequent evo¬ 
lution histories. The flattening and high dispersion of the AMR in 
the solar neighbourhood is often ascribed to effects of radial mi¬ 
gration in the disc due to non-axisymmetric perturba tions of the 
gravitational potential (e.g. [Sellwood & Binnevll2002h . this region 
containing stars born in a wide range of galactocentric distances. 
The increasing AMR in the outer disc of M31, as well as the com¬ 
plex and differing SFH behaviours observed in the three fields, 
suggests that migration in this galaxy, or at least this part of the 
galaxy, has been less efficient than in the solar neighbourhood. For 
example, it is likely that the old, most metal-poor component of 
the Warp would have had sufficient time to contaminate the 23 kpc 
field had migration been significant. The different migration effi¬ 
ciencies may be related to the different dynamical properties of the 
M31 and Milky Way discs: IVera-Ciro et al.l(l2014l) have shown that 
the efficiency of radial migration is a strong function of the verti¬ 
cal velocity dispersion, the stars spending more time near the disc 
plane being more affe cted by the non-axisy mmetric perturbations. 
The measurements of IPorman et al.l 1 I 2 OI 5 I) in the disc of M31 in¬ 
dicate a velocity dispersion at a give n age that is more th an twice 
as large as in the Milky Way (see also lCollins et al.l ( 1201 ih l. poten¬ 
tially providing an explanation for the lack of migration. The old 
age of the M31 outer disc provides further indirect support for this 


idea, since disc galaxies observed at 2 are often highly turbulent 
(e.g. lFdrster Schreiber et alj2009l) . 

In summary, our analysis of deep CMDs that reach back to the 
oldest MSTOs have enabled the first precision constraints on the 
radial variation of the star formation history in the M31 outer disc. 
Roughly half of the stellar mass in our fields was formed before 
z ~ 1, indicating that the extended M31 disc was in place and ac¬ 
tively forming stars at high redshift. Stellar feedback was not strong 
enough to significantly delay the formation of the outer disc. Our 
results also suggest that radial migration has not played a major role 
in rearranging the stellar populations at large radii in M31. Notably, 
a well-defined AMR is seen in all fields, in marked contrast the flat 
AMR seen in the Milky Way at similar galactocentric distances. 
Our findings provide new insight on processes relevant for mod¬ 
elling the evolution of disc galaxies, issues which we will explore 
in more detail in future work. 
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